Meson Cloud, QCD Nucleon Correlator and Chiral 

Perturbation Theory 

Leonard S. Kisslinger 
Department of Physics, 
Carnegie Mellon University, Pittsburgh, PA 15213 

February 1, 2008 

Abstract 

We introduce a model for nucleons with correlators that include the prop- 
agation of pseudo Goldstone bosons as well as the standard correlators for 
nucleons used in QCD sum rules. Prom the comparison with experimental nu- 
cleon magnetic dipole moments we estimate the probability of a pion cloud in 
a nucleon. The masses obtained from the sum rules are almost independent of 
the cloud contributions. We show that the model is consistent with chiral per- 
turbation theory, and solves a long-standing problem on obtaining the leading 
non-analytic correction to the nucleon mass in QCD sum rules. 
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1 Introduction 



The concept of a meson cloud is very old in nuclear physics: the existence of the 
pion exchange interaction leads to a pion cloud [[L| which could affect the properties 
of nucleons. Although a main current goal of hadronic physics is to understand the 
properties of hadrons in terms of QCD, in recent years the study of deep inelastic 
scattering (DIS) has led to a great deal of interest in the nucleon's meson cloud. 
The NMC/CERN experiments^ gave evidence for the violation of the Gottfried sum 
rule and the recent Fermilab/E866[|3], |J Drell-Yan measurements show that for small 
Bjorken x the ratio d(x)/u(x), down to up sea quark distributions, is considerably 
larger than 1.0, while both perturbativeU and nonperturbative QCD calculations || 
show that d(x) / u{x) for small x differs from unity only by about 1 % isospin violations. 
It is likely that theoretical treatments with both QCD and Goldstone boson (GB) 
fields are most efficient in treating these sea quark problems as well as problems 
with the spin and strangeness content of the nucleon. The challenge is to do this 
consistently. 

There have been many theoretical calculations of DIS using the concept of a meson 
cloud based on the one-pion interaction of nucleons 0] (see Ref.[§ for a review). It is 
obvious that if one only considers the 7r + -neutron system as an additional component 
of the proton, that there are more d than u sea quarks, but the 7r~A components 
tend to cancel this effect. A related model is the chiral quark model which has 
been used for flavor and spin properties of nucleons|9|]. As discussed in Ref.||, al- 
though these two models both introduce Goldstone bosons, in the chiral quark model 
these are fields internal to the confinement region, while the meson cloud model de- 
scribes long-distance effects outside of the region of confined quarks. In lattice gauge 
calculations |TtJ of the flavor-singlet axial coupling constant for the nucleon's quark 
spin content gluonic effects were included by disconnected quark loops, which are 
interpreted as representing Goldstone bosons. In this method valence quark and GB 
effects are distinguished. 

It was observed that the QCD sum rules for the nucleon mass are inconsistent 
with chiral perturbation theory [fTT|1 . First, the quark condensates terms on the QCD 
side of the sum rules contain quark condensate terms which have been shown Hl2[ 
to introduce m^ln(-m^) terms, while the chiral perturpation theory expansion for 
the nucleon mass[|l3|] does not have such terms. Second, the leading non- analytic 
correction (LNAC) to the nucleon mass is a m^/fj; term, which is not found in QCD 
sum rules IP- By an analysis of the phenomenological correlator, given by a term with 



a pole at the nucleon mass and the continuum, it was shown WM that the chiral logs 
cancel; however, the LNAC term is not present so the QCD sum rules for the nucleon 
mass are not consistent with chiral perturbation theory. Virtual pions have also been 
included in QCD sum rule calculations of nucleon masses in nuclear matter [|15|, I6|. A 
chiral quark model was used to study pion-nucleon weak and strong coupling within 
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the QCD sum rule method 

The model introduced in the present work uses the method of QCD Sum Rules |18| 



in which hadronic systems are represented by local complex field operators, usually 
called currents, so that two-point functions can be used for defining correlators from 
which one can determine hadronic masses. Nonperturbative effects are treated via 
operator product expansions (O.P.E.) of the quark and gluon propagators using vac- 
uum condensates, whose values are determined by fits to experiment as well as lattice 
gauge calculations. In this method short-distance effects are treated perturbatively 
and the condensates account for the medium-range and long-range QCD effects within 
the region of confinement. 

The model introduced in the present work treats the nucleon correlator with a 
nucleon current that has a quark/gluon component and a component also involving 
a Goldstone boson field. The GB fields are included in the composite-field nucleon 
current, with no coupling to either hadrons or quarks and thereby represent the long- 
distance meson cloud effects. Thus the model differs from both the usual meson cloud 
and the chiral quark models. All internal meson exchange effects are assumed to be 
included in the condensates. For consistency, the coupling between the two currents is 
also not included, so the correlator has a standard QCD sum rule part and a QCD / GB 
part, with a GB propagator. There is a clear relationship between the present model 
of nucleons with a meson cloud and the usual QCD sum rule treatment. In the usual 
treatment states with a nucleon and a pion are in the continuum. In the present 
model the current explicitly contains pionic fields, so that the resulting continuum 
is modified. By doing this we are able to treat physical effects where explicit pion 
degrees of freedom are important, as in sea-quark distributions. 

The goals of the present work are to estimate the probability of the meson cloud 
component by studying the properties of the nucleon, and to study the problem of 
consistency with chiral perturbation theory. This is similar to the problem of finding 
the correlator for hybrid hadrons JOI, where the hybrid components of the nucleon 
and lowest hybrid baryon were estimated from the hadronic masses. In Sec. 2 we 
introduce the correlator using a standard form|2(J for the nucleon with no GB and 
derive the sum rules for the nucleon masses including a correlator with a GB. We 
find very small effects from the meson cloud correlator terms. In Sec. 3 we derive 
the sum rules for the neutron and proton magnetic dipole moments using a two-point 
correlator in an external electromagnetic field From the sum rules found with 
this method compared to QCD sum rule calculations of Ref. |21| without the cloud, we 
conclude that the GB component has roughly the same probability as the cloudless 
component. It should be noted that there is recent work on including instanton effects 
on the quark propagator that also modify the sum rules for the nucleon moments p2 



In Sec. 4 we show that the model is consistent with chiral perturbation theory. The 
chiral logs cancel and the LNAC term appears as in hadronic models. The known 
value of the LNAC term from chiral perturbation theory can help determine the 
parameters of the model. 
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2 Correlator with Goldstone Boson Fields: Nu- 
cleoli Masses 



In this subsection we present our model of the pion cloud within the QCD sum 
rule method, which means that we introduce a correlator with explicit Goldstone 
boson fields. We apply this model correlator, a two-point function, to the estimate of 
the masses of nucleons. The composite field operator for the proton has the form 



■nix) = Cl 7f° + czrf'* + c 3 ri P > K , (1) 



with constants Cj giving the meson cloud amplitudes. In Eq.(|I]) r] p '° is a standard pp] 
proton field operator used in sum rules: 

rjP>°(x) = e abc [u a {x) T C^u b {x)]j%d c {x), 
< 0\r] p '°(x)\proton > = X p v(x), (2) 

where C is the charge conjugation operator, the u(x) and d(x) are quark fields labelled 
by color, X p is a structure parameter and v(x) is a Dirac spinor. The neutron operator, 
rf 1 ' , is obtained from r] p '° by interchange of u and d quarks. The composite field 
operator r) p ' w used in our present model is 



= -i_[^ Q 0+(x) 7 a 7V'° + <9 Q C(x) 7 Q 7V 
< $\ri(x) p '*\proton > = X' p v(x), (3) 

where 4>tt(x) is the pion field and is a D=l scale factor. From our experience 
with hybrid baryons [fl9[| and from the p-wave coupling of the pion one expects that 
X p 2 <C Xp. This simplifies the phenomenological form for the correlator and enables 
us to obtain sufficient sum rules to estimate the pion cloud part of the current. In the 
present work we do not consider the kaon cloud term which would enter through rj p,K , 
defined similarly to Eq. (|3]) but with K fields and strange baryon composite operators. 
For estimates of strangeness in nucleons, which is a natural application of our model, 
this term must be included. 

Our resulting model correlator for the proton with a pion cloud is 

IP (p) = i J d 4 e txp < 0\T[ri(x)r](0)]\0 > 

= cli J d A e lxp < 0\T[r] p '° {x)ff '°{0)]\0 > +(1 - c\)i J d 4 e lx ' p < 0|T[7^(x)7f >"(0)]|0 > 
= clU {p '°\p) + (l-cl)U^\p). (4) 
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An important aspect of the model is that the coupling of the n p ' n and n p, ° currents are 
not included. Only the long-range effects of the pion are included by the pion cloud 
contribution, and no pion-quark interactions are included in order to avoid the double 
counting of interactions given by the condensates. In a chiral quark model the pion- 
quark coupling would give both the coupling of the two currents and internal pionic 
forces. The present model assumes such internal interactions are given by QCD. 
Only external fields interact with the pion, giving the physics of the meson cloud. 
The neutron correlator can be obtained from Eq.(|4]) with obvious isospin changes. 
The correlator Jv 3> '°^(p) for the propagation of a proton without the meson cloud is 
given in Ref.|21|] and many other references. 

Next we discuss the correlator with the pion cloud, U^{p). From Eqs.(|.|) one 
finds for U^\x) 

n^C*) = ^ abc e b,a,c, Gi,{x) lal ^'{x) lulp 

Tr[Sf(x)YC(S b u b '(x)) T CY] 

+2[n° — > 7r + and u <-> d] + [4 — quark terms], (5) 
where in the limit of zero pion mass 

W = i(^-M- (6) 

In momentum space the pion cloud correlator is 

TJMM 9 f d ' k (P-m P '°(k)(P-k) 

[p) ~ J (2tt) 4 (p-ky 

= pU^\p) odd + U^\p) even , (7) 
with p— r ) c p a . We find after the Borel transform from p 2 to the Borel mass variable, 

2 
B 



M 2 



B ! even 



(27r) 4 n^(M 2 ) odd = _ 2§ [12M^°E 4 - 5bM B E 2 - 40a 2 M%E X 



+20mia A M£E ), (8) 

where the parameters for the quark condensate, gluon condensate and mixed con- 
densate, respectively, are a = -(27r) 2 < qq >, b = < qg 2 G 2 q > and m 2 a = (27r) 2 < 
qga ■ Gq >. For these parameters we take a = .55 GeV 3 , b = .47 GeV 4 and m 2 = 
0.8 GeV 2 . The scale factor A 2 has a value in the range of M 2 B /g nN (see Sec. 4). 
The functions E n are l-exp(-w)X^=o wfc A!) with w = s /M 2 B , s a being the threshold 
parameter. In comparison with Eq.(||) for the meson cloud terms in the correlator, 
the usual(p,0) odd correlator is 

(2vr) 4 n( p '°) (M 2 B ) odd = ^(-4M B E 2 + bM 2 E + lQa 2 /3-4m 2 o a 2 /(3M 2 B ). (9) 
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The dependence of the correlators on the anamoulous dimensions are not shown in 
Eqs.(|,^|), since the result is not essentially changed if this is neglected. In fact it 
is obvious from Eqs.(||,fJ) that if A p 2 << A 2 , as expected, then the sum rules for the 
masses are independent of the pion cloud within the accuracy of the method. 

3 Correlator with Goldstone Boson: Nucleon Mag- 
netic Dipole Moments 



For the calculation of the magnetic dipole moments of the proton and neutron we 
use the external field formalism, which is one way to avoid the problem of carrying 
out an operator product expansion in the momentum transfer variable. The coupling 
of a current J T to a nucleon is given by a two-point correlator in the external J T 



current 21 



n r (p) = z J d 4 e tx p < 0|T[t7(x)t7(0)]|0 > jr . (10) 

For the proton in an electromagnetic field, represented by the field tensor F M1/ the 
correlator is written as 

= e{c\W£ + (1 - cl)U^)F^, (11) 

with the two terms representing the correlator in the electromagnetic field without 
and with the meson cloud. Since at low momentum transfer the coupling of the 
photon to the pion does not contribute, the meson cloud correlator in the presence of 
the field in our model is given by 

^ [p) ~ J (2tt) 4 (p-ky 1 ] 

There are three covariants, and we write the correlator with and without meson clouds 

as 

IW*;) = [a^,fc]nf dd (A ; ) + z(^7,-^7M)nL n (A;) + a^nf(A;). (13) 

For the two-point effective field treatment at low momentum momentum transfer 
n r (:r) can be evaluated using the operator product expansion (O.P.E.), since the 
variable x is at short distance from the origin for large p M . This is done by an O.P.E. 
of the quark propagator in the presence of the the J r current 

S T q (x) = <0|T[g(x)g(0)]|0> jr . (14) 
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There are three susceptibilities for the induced condensates: x-, K an d £, which are 
defined in Ref . For example, the quark condensate magnetic susceptibility, X , is 
defined as 

(27r) 2 < qa^q > F = -e qX aF^ v . (15) 

The QCD evaluation of the correlator in an electromagnetic field without a meson 
cloud gives pi| 



(27r) 4 n2£ F (A:) = e u k 2 ln(-k 2 )-e uX a 2 /(3k 2 ) + C a /(3k 4 ) 

(2ir) 4 W e ^(k) = a[(e u + e d /2)/k 2 + (e dX /3)(ln(-k 2 )-b/(24k% (16) 

with C a = {a 2 /2)[e d + 2e d /3 - e u (n - 2f)] - e u xa 2 m 2 j8. From Eqs.(g[RjD we obtain 
the correlator for the nucleon with a pion cloud in an external electromagnetic field, 
giving 

(2*)*I&(p) = l -^^[-euP 6 + ^e uX aY/3 + 32C a ] 
WU^(k) = - ^t^_ xa[3 K 2 ^-5e d b\. (17) 

For a rough estimate of the effect of the meson cloud, which is the goal of the present 
paper, we follow the prescriptions of Ref.ETfl: For the odd term in the correlator 



multiply the p term by and subtract e u times the n term. Neglecting anamolous 
dimensions one finds 

—A 2 a 2 
-^e- M ^ M B(e df i p -e u fi n )/M 2 + single poles = (e 2 - e*)[^j + A d ], (18) 

where A^ = (27t) 4 Aat and the GB contribution is 

A c , = \ ~ A f 27 ^ 2 [8a 2 M 2 B E - 6M|E 3 - f X a 2 E, + |a 2 E + ^a 2 m 2 ].(19) 

Noting that with a value of x a = "4 GeV, which we obtain with our three-point 
treatment with nonlocal condensates^], a value of c 2 about 0.5 reduces the neutron 
and proton magnetic dipole moments by about ten percent, improving agreement 
with experiment. From this we conclude that the component of the correlator with a 
pion cloud is roughly equal to that without. As we shall see in the next section, the 
known LNAC for the nucleon mass will help determine the scaling parameter, A w . 

4 Consistency of Model With Chiral Perturbation 
Theory 



7 



In Ref. [0 it was shown that from the dependence of the quark condensates 
on the current quark masses and the Gell-Mann-Oakes-Renner relationp4j that the 



quark condensate has a m^n(-m^) dependence on the pion mass. From this one can 
see from Eq. that the standard QCD sum rule contains chiral logs that are not 
consistent with chiral perturbation theory Moreover, the LNAC term, which is 



proportional to m^/f^ |13] does not occur in the QCD sum rule form for the nucleon 
mass, in this section we show that both of these problems are solved in our model. 

4.1 Chiral logarithms 



In this subsection we show that by using the results of Ref. [[L4j] our present 
model does not contain inconsistent chiral logarithms. We study the diagrams giving 
the LNAC in the next subsection. For simplicity let us take c\ = 0.5. Using the 
observation from hybrid baryons[|iTf| that A p 2 << A 2 , as discussed in Sec. 2, the 
phenonenological side of the correlator, given by a dispersion relation has the form of 
a pole term and a continuum term, 

2W( phen \p) = -\ 2 / 2 _ J 2 +W^ nt \p), (20) 

V p 

which is the same as the standard proton correlator except for a factor of two. As we 
saw in Sec. 2, the microscopic QCD side of the correlator has the form 

2w^ d \ P ) = p(n^ ) QOT (p) odd + n^ QOT (p) odd ) + n^°)(p) et)en . (21) 

Using the standard proton current we recall that 

b 

18p~ 2 



U^ QCD (p) even = < q q> (p 2 ln(-p 2 ) -—) + ■■■ (22) 



n {pmcD (p) odd = -H- P 2 )(^ + b -) + ---, 

where by the ellipsis we mean four-quark condensates + higher dimensional terms. 



Making use of the m n expansion of the quark condensate from Ref. [|T^] : 

< qq > = < qq > a [1 — cm 2 Jn(— m 2 )] (23) 

where c = 32 n 2 / 2 anc ^ the symbol Q Q we mean the quantity Q in the limit m 2 — > 0, 
one finds 

U^ CD (p) even = [l-cmlln(-ml)}U^ CD (p) even + 0(ml) (24) 

n(pmcn {p)oM = [i + (ml)]rtr )QCD (p)oM, 

as was observed in Ref. [T^] . For the pion cloud part of the correlator we observe that 

n.M(x) = G^(x)7 a n>°y + m^G^(x)7 Q n>°y, (25) 
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where G a p is obtained from G Q/3 of Eq.(^) by an expansion of the pion propagator in 
powers of m 2 , with the first term given by Eq(^|). From Eq.(|25| we observe that 

ILM(x) = n<?ri(x)+0(ml), (26) 

and does not contain any m 2 ln(—m 2 ) chiral log terms. 

Therefore we conclude that in our model of the nucleon with a meson cloud 

W [Phen) = [1 _ cm 2/ n (_ m 2^ n p(Pft e n) + . . . ^ 

and 

W (QCD) = [ 1 _ CTO 2 Zn( _ m 2pp(QCi3) + .... (2g) 

and our model the nucleon with pion cloud does not contain spurious chiral loga- 
rithms. 

4.2 Leading Non-Analytic Correction 

The leading non-analytic correction (LNAC) from chiral symmetry breaking can 
be found in chiral perturbation theory using the method described in Refs. 13 



with the Feynman-like diagrams of chiral perturbation theory given in Ref.[26]. The 



main idea is that if one has a broken symmetry, so that the Hamiltonian has the form 

H = H + mH u (29) 

with m being the small parameter of the symmetry breaking, the mass of a hadron 
can be found from the relation 



d 2 M 2 
dm 2 



Lmv^o+^o.p^ / dSe iqx < p|^(x o )[H / (^),^/(0)]b >, (30) 



including only connected diagrams, with Hj = / d 3 x7ii. The chiral symmetry break- 
ing in QCD is 

mHi = J d 3 xm u uu + niddd 

= m f d 3 xH m . (31) 



In Ref |13| it is shown that the LNAC for the nucleon mass is given by the 
process illustrated in Fig. 1. The corresponding correction to the nucleon mass, 
A = Mat — was shown to be 

AM N = -Ml§ + ... = -i 5 MeV + ---, (32) 



71 



Or 



\ 



N 



Fig. 1 Diagram giving LNAC to nucleon mass 



which makes use of the result 

2 

<7i a \H m \7c b > = 5 ab ^ + ---. (33) 

m 

The LNAC in the present model of QCD sum rules, after the chiral logs have been 
eliminated, are given by processes involving the pion cloud correlator. The lowest 
dimensional term is illustrated by the the diagram shown in Fig. 2 In the evaluation 

71 

Or ~ P 




Fig. 2 QCD process giving LNAC to nucleon mass 



of the QCD sum rule diagram the 7r-N T-matrix evident in Fig. 1 is replaced by 
the traces over the quark propagators. The LNAC to the nucleon mass is found by 
equating the phenomenological expression for the Borel transformed correlator 



M 2 

U p(Phen)^ = p2 exp {__JL} +continuum ^ (34) 



with j3 2 = (27r) 4 Ap/4 (see Eq.@, to the QCD expression with chiral symmetry break- 
ing. Making use of Eq(|33"D, the Goldberger-Trieman relation (g 7r iv=g J 4 M N /f T ), Eq.(|j) 
and taking the nucleon and Borel masses M p , Mb — 1 Gev, we find for the lowest- 
dimensional LNAC term 



AM p 



em\g\ 



-20MeV. (35) 
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For the result of a LNAC of about 20 MeV, which is consistent with the chiral 
perturbations theory result, we have used a value of (3 2 from Ref. [21], where it is 



pointed out that there is about a 50% uncertainty in the value, and have used = 
Mp/g^Tv for the scale factor of the pion cloud model. Note that the gluon condensate 
process also contributes to the LNAC, with an opposite sign, cancelling about 10-20% 
of the term shown in Eq. (|35|). 

Therefore, we have shown that the longstanding puzzle of the absence of the LNAC 
m^ term in the QCD sum rule treatment of the nucleon mass can be explained in our 
model, with a correlator composed about equally of a conventional quark correlator 
and a meson cloud correlator. The known value of the LNAC for the nucleon mass 
can help determine the parameters of the model. 



5 Conclusions 

We have introduced a model for a nucleon in which Goldstone bosons are included in 
a QCD treatment. By only including the boson propagator and not interior quark- 
boson interactions we hope to account for long-distance effects which are difficult to 
calculate in a pure QCD treatment, and yet avoid introducing processes which are 
aleady treated by the nonperturbative QCD methods of QCD sum rules. With this 
model the nucleon mass sum rules are essentially not changed by the meson cloud 
terms, but from the treatment of magnetic dipole moments we conclude that within 
our model the correlators with and without the pion cloud are about equal. We 
have also shown that this model is consistent with chiral perturbation theory, and 
that the leading nonanalytic term found in chiral perturbation theory, but missing 
in previous treatments of the nucleon mass with QCD sum rules, is present in our 
model. The magnetic dipole moments and the LNAC for the nucleon mass can help 
determine the parameters of the model. Our future research with this model will 
include studies such as the spin and strangeness content of the proton, as well as the 
sea-quark distributions in the nucleon. 

The author would like to thank Montaga Aw, Andrew Harey, Ling-fong Li and 
Jen-Chieh Peng for helpful discussions. This work was supported in part by NSF 
Grant PHY-9722143. 
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